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O ■ Abstract 

O I Large order asymptotic behaviour of renormalization constants in the minimal subtraction scheme for 

the (j)'^ (4— e) theory is discussed. Well-known results of the asymptotic 4— e expansion of critical indices 
are shown to be far from the large order asymptotic value. A convergent series for the model (/)^ (4 — e) is 
then considered. Radius of convergence of the series for Green functions and for renormalisation group 
p • functions is studied. The results of the convergent expansion of critical indices in the 4 — e scheme 

are revalued using the knowledge of large order asymptotics. Specific features of this procedure are 
discussed. 
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H 1 Introduction 



Calculation of critical indices is usually based on a certain asymptotic expansion. To obtain reliable 
results a resummation procedure is necessary. To this end the Borel-Leroy transform in fixed dimension 

^ or in e expansion ||^, ^, the simple Pade-Borel method ^ as well as self-similar exponential 
approximants ||^ have been used. But divergent series can produce an arbitrary result. Thus some 
additional information about the series is needed. 

When a convergent series is considered, additional information can also improve the result. The 
knowledge about location and character of singularity of the function investigated or about the large 
order asymptotic behaviour of the expansion may be used to accelerate convergence of the series. In this 
report we will discuss a resummation of the 4 — e expansion of critical indices for the (p'^ theory with 
help of the large order asymptotics. 
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2 Large order asymptotics of renormalization constants 



Large order asymptotic behaviour of the 0(n) symmetric (j)^ theory is well known. The N-th order term 
of the standard expansion for an arbitrary function F in the theory with the renormalized action 



Sr{4>,9) 



i2\2 



behaves at large N as (O, § m 



(1) 



(2) 



The notation F^^^ for the A^-th order coefficient of the expansion of the function F will be used hence- 
forth. Working within the minimal subtraction (MS) scheme in the 4 — e expansion we have corrected 
the result (^. Namely, analysing the limits N ^ oo and e ^ in more detail, we have obtained a 
more accurate estimate for the amplitudes Ci? of the renormalization constants [[To[]. According to the 
results of Ref. the asymptotic expression @ for the /3 function coincides practically with the exact 
result for = 3, but the asymptotics of [^] for = 4, 5 are larger than the exact value. Contrary to 
these conclusions, we have shown that the asymptotics for the renormalization constants (and for the 



critical indices) are much smaller than the exact ones Jl For example, [Zg ]asymp ~ •01[Z, 
n = 1, where [Zi] is the residue of Z,; regarded a function of complex e. We have predicted that only 
starting from 10 - 15-th order the perturbation expansion could be near to the asymptotic value. Due to 
this fact we can state that the Borel transform of the 5 known terms of the e expansion has no theoretical 
ground. 

Nevertheless, we propose extrapolation expressions for the unknown terms in the expansions of 
renormalization constants Z in the following form [ ]T0| ] : 
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These expressions contain an additional correction of the type normalized by direct comparison of 
the asymptotic value (^ for = 5 with the exact one. In (||), 
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^'(1) - ln(7r) 



^'(1) - lii(7r) 
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nr(n/2) \ 6 

jj^^ ^ 2(''~''")/23W2^-(5+n)/2g-5/2g-4(n+2)/3-7(n+8)/18-fl/12+(n+8)(1'(l)+2-/n(7r))/6 

Here, is the logarithmic derivative of the F function. Using similar expressions normalized with help 
of the 4-th exact term of the renormalization constant expansions, we could find the 5-th order term and 
estimate the accuracy of such calculation as at least 80%. Besides, we hope that the expressions (^) can 
give us a reliable estimate for more than 5 orders of the expansions of Z's. 
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3 Convergent expansion for critical exponents 



We will use the knowledge of large order asymptotics to improve initially convergent series for critical 



exponents. In Ref. [12] a new approach was suggested to calculate critical indices. It was based 



on the modified expansion for the 0(n)-symmetric 0^ model [13] which leads to convergent series. 
The standard action S = Si + S2 with Si = f dx{d(j))'^ /2, S2 = g f dx(p^/4l was cast in the form 
S = Sq + Sj, where Sq = Si + aSl and 5/ = C(5'2 — aSl). Here, a is an arbitrary constant and ^ is a 
new expansion parameter. The model coincides with the initial one at the 'physical' value of Q: C,p}i = 1. 



It was stated in Ref. [13] that the C, expansion is convergent for < = 1» when 



Q- ^ g47i-2 — ^min- (7) 

In the framework of the C, expansion the renormalized 2k point Green function can be written as an 
integral with respect to an additional variable a as: 



J {l + 2ia^a{l-Q))^ ^ 

X expf-lw C9^^\_ 

2 4!(l + 2iav/^(W))2y 

Here, g is the renormalized charge, /i is the renormalization mass, Zi = Zi{g(/{1 + 2ia\J a(l — C) )^), 
and Zii^g) are the renormalization constants in the usual perturbation theory. This representation allows 
to construct the Feynman graphs in the usual manner. Due to the integration with respect to a every 



diagram acquires an additional factor [ ]13[ ] : 

U-[a) = l/V^ j daexp{-a'^)/{l + 2ia^/^y 

dtt^^^exp{-t - at^)/{j - 1)! (9) 



The basic renormalisation group (RG) equation for the Green functions (g) was derived in Ref. [ ]12| ] in 
the form 

T^RGG^k = + pdg + (32dl + ... + 7]Gffc = 0, (10) 

where 7 and all (i functions depend on the parameters g, C,, a and k. Using the MS scheme one can write 
for the RG functions in D = 4 — e dimensional space (n = 1, k = 1) the expressions 

2 « lAfyW 



7 = - Eg' 1^2^+1 ("(1 - ^))' (11) 



= -eg + l3{g)g (12) 

It was shown that Eq. (pX|) governs the large-scale asymptotic behavior of the model. Critical exponents 
are related to the anomalous dimensions in the usual way [^: r? = 'y{g*)/k, and l/u = 2 — ^r{g*) + krj 
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[ 7r is determined by equation ( |10| ) for the Green functions with the insertion of the composite operator 
(f)^ ]. Here gf* is the fixed point determined by the usual equation /3(5*) = 0, or 



(13) 



In Ref. [12] it was demonstrated that it is possible to solve Eq. (13) iteratively to calculate g^, in the form 
of a double expansion in e and Q. As a result, the exponents rj and v were calculated forn = 1, A; = 1, 
a = 0.134 at the physical values Q = e = 1 (see the table). The columns correspond to expansion order 
taken into account. Exponents marked by e are the results of the usual e expansion, which are quoted for 
comparison. 



Exponent 


1 


2 


3 


4 


5 


V 





0.02553 


0.04342 


0.039745 


0.039741 


Ve 





0.01852 


0.03721 


0.02888 


0.05454 


u 


0.621 


0.671 


0.663 


0.674 


0.651 




0.583 


0.627 


0.607 


0.678 


0.461 



Estimates based on the Borel transform of e expansion are -q = 0.0360 ± 0.0050 r] = 0.035 ± 
0.002 [|], V = 0.6290±0.0025 [|], v = 0.628±0.001 [§. A typical lattice resuU is = 0.6305±0.0015 
[^. Thus, in spite of convergence of the C expansions, the accuracy of the results of [ p^ ] is lower than 
in case of estimates based on Borel transforms. This is why we will improve the convergent expansion 
using the large order asymptotics. 



4 Large order asymptotics for the convergent perturbation 
expansion 

In Ref. [ [I3| ] it was stated that instanton analysis could not be used for the investigation of the convergent 
(^-expansion of Green functions. Using more sophisticated examination, the radius of convergence was 
estimated as = 1. 

Contrary to this, we confirm the adequacy of instanton approach here. Indeed, let us determine the 
A^-th order of C, expansion by 

.(TV) _ I dC 



^2k - J ^lvTT^2fc, 

where G2k is given by (|8|). A subsequent steepest descent approach in the variables a, ( and (f> leads to 
an instanton and stationary point (for n = 1) 



where 

4^/3 



y l + (x-xo)2/y2' 

with arbitrary xq, y, and 



6V^ 'N'r l-5/(647r2a) 
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AT U I ^ , SeSavr^ + 35 , 1 . 

"° = V C^V^ + 96vraVg + 

The instanton analysis results in the following behavior for the A^-th order of the (^-expansion of the 2k 
point Green function: 

G(f -iV^Cfe-v/I (14) 

with some a. 

Contrary to [ [D] ] we have obtained a radius of convergence which tends to infinity as a ^ dmin- 
We consider the unrenormalized Green function in D = 4 dimensional space only, but this is not essen- 
tial for the treatment of the radius of convergence. The same result (0) can be obtained directly by the 
method proposed in 1 13]. In such a case the instanton 0c ~ N'^^^. Thus, the instanton does exist. 



For usual divergent series of e expansions the large order asymptotics of Green functions characterise 
unambiguously the asymptotics of RG functions. In the case of convergent series the situation is more 
difficult. Namely, the denominators of the expressions ([iT]), ([T2|): Ui{a{l — Q)) {i = \, 3) have the radius 
of convergence equal to 1. Therefore, the large order asymptotics of the C, expansion of these functions 

- 0) ^ ivV2-ie-^+i/(8a)^i-Yp(.) 
have a non-trivial influence on the asymptotic behaviour of /3 and 7 functions. 



5 Character of singularity of the convergent expansion for 
critical indices 



In addition to Q problem we have some difficulties with double C> e expansion of indices as we try to 
find the fixed point of (|l3]). The iteration solution fif*(e, C) of the equation (|T3]) has a new singularity 
at the point bounded by the nearest to zero extremum of the j3 in the complex plane g. Namely, the 
point gm given by 



dgPign 







(15) 



results in the singularity g^ ~ — e^. 

To calculate em the second equation (|l^) has to be solved in a form of a C expansion. Unfortunately, 
our calculation shows that the five known terms of the expansion of /3 are not sufficient for a reliable 
determination of gm- If the extrapolation expressions (js]) are used within the Pade approximant approach, 
then there is a whole set of conjugate points g^- It is very difficult to find a suitable 'mapping' to exclude 
these singularities. Moreover, the corresponding minimal value of turns out to be very small (~ .02) 
that decreases the radius of convergence in e of our double expansion. Thus the method to be used is to 
consider the variables g, ( instead of e, C- A similar approach leads to a good result in the Kraichnan 



model, where a convergent series is dealt with as well [14]. 



Calculating g^ directly from the equation ( |13| ) for e 
expansion for critical indices (for n = 1, (ph = 1): 



1 we were able to obtain five orders of the C 



Exponent 
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2 


3 


4 
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V 





.021 


0.023 


0.027 


0.029 


5*/ (16^') 


0.96 


0.86 


0.65 


0.61 


0.57 


V 


0.612 


0.650 


0.638 


0.641 


0.640 


5*/(16^') 


1.38 


1.17 


0.84 


0.75 


0.68 
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The exponent r/ was calculated for the value a = .14, the exponent u for a = .17 which lead to the 
best rate of series convergence. Note that obtained ensures the convergence of the C, series in both 
cases. However, using the extrapolation expressions (|3|) for unknown terms of the expansions of the 
renormalisation constants Z, we obtained the results shown in Fig. 1 demonstrating the failure of the 
usual extrapolation procedure in the calculation of critical indices. 

To improve the convergence of the ( series we can investigate the large order asymptotics of the 
functions (pTj), (|I2|). Denominators of these expressions have essential singularities of the type Ui{a{l — 



(), U3{a{l — C) that must be extracted. It is convenient to introduce SZ^^ = [Z^^] — [Z^'^]asymp- 
Let us rewrite the asymptotic expression (Q) of the renormalization constant Z^p in the equivalent form 

Here, the constant can be found in a way similar to in @ by comparison of [Z^^]asymp with the 

(5) 

exact expression [Z\ ]. 



Then substituting [Z^^]asyrnp into (11) and re-expanding it in ( one can write 7 in the form 



^ (16vr2V Ui{a{l - ()) C/i(a(l - C)) 

^ ( - p^^'^^ + (^''■^' + - (B-f^' + ^^)^^) • ^''^ 

Here, the omitted for brevity argument of the Ui functions is a(l — C) + 5C/647r^. Thus, we have 
investigated not only the location of Cc and Q, but also the character of the singularities of the C expansion 
of the 7 function. 

In analogy with ( [T^ ) we obtain for the /3 function 

. ^ . {gCy U2H3Hl-0){SZ^g'^-26zl^^) c, 

jr{\^Q^^V ^3Ki-C)) +25/2[/3(a(i-c)) 

( ^90' g,^g ^ 7i?^6!(3 + Cg)U'r - 7Mt4!(1 + Cg)uA - 7 (17) 



I (647r2)3^-^« (647r2)2"-^" ^3'^' (647r2 

The expressions ([I^, ( |l7| ) are physically meaningful at = 1. Solving numerically Eq. ( [I3| ) with fi 
from (17) and substituting the value g^ = 0.382 obtained in 7 (p^) we obtain rj = 0.0236. Taking into 
account the e expansion up to fourth order one obtains similarly 7^4 = 0.0241. 

An analogous procedure for the index v leads to = 0.580, (1/4 = 0.624). These results are the best 
we can obtain from the e expansion in the convergent scheme using all available information about the 
large order behaviour. 

It is worthwhile noting that the accuracy of our results is similar to the accuracy of the approxima- 
tion (|3|) determined by the rate with which asymptotics of renormalization constants tend to the exact 
value. Thus, we conclude that the accuracy of the e expansion resummation with help of the Borel 
transformation is lower than commonly quoted. 
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Figure 1 : The value of the index 77 and the index z/ as a function of perturbation order number taken into 
account. 



8 



